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RADIOMETRIC  TECHNIQUES  EOR  OBSERyiNG 
THE  ATMOSPHERE  FROM  AIRCRAFT" 


William  L.  Smith  and  Warren  J.  Jacob 
National  Environmental  Satellite  Service 
National  Oceanic  and  Atmospheric  Administration 
Washington,  D.  C. 


ABSTRACT.   Radiometric  observations  have  been  made  from 
aircraft  with  spacecraft  prototype  instruments  to  test 
satellite  remote  sensing  techniques.   At  the  same  time 
these  observations  have  been  used  to  study  the  problem 
of  remote  sensing  from  aircraft  because  of  its  own 
particular  value  in  providing  data  in  the  support  of 
sub-synoptic  scale  meteorological  experiments  such  as 
the  forthcoming  GARP  Atlantic  Tropical  Experiment. 

This  paper  describes  aircraft  radiometric  methods  of 
obtaining  clear  column  radiances,  vertical  temperature 
profiles,  total  precipitable  water,  and  cloud  heights 
and  amounts.   Questions  regarding  vertical  resolution 
and  accuracy  specification  as  a  function  of  aircraft 
altitudes  are  answered.   Results  obtained  using  obser- 
vations made  during  June  1970  with  prototype  versions  of 
the  Nimbus-E  ITPR  aboard  the  NASA  CV-990  are  presented 
and  discussed.   It  is  concluded  that  multi-spectral 
radiometers  aboard  an  aircraft  can  be  effective 
meteorological  observing  tools . 


1.   INTRODUCTION 

In  this  paper  we  present  methods  of  making  meteorological  inferences  using 
multi-spectral  infrared  radiance  observations  made  from  an  aircraft.   In 

particular,  methods  are  given  for  obtaining:   (l)  clear  column  radiances, 
(2)  the  vertical  temperature  profile  from  aircraft  altitude  to  the  earth's 
surface,  (3)  the  total  precipitable  water  beneath  the  aircraft,  and  (4)  the 
height  and  amount  of  any  clouds  below  the  aircraft  flight  level.   The 
techniques  are  tested  using  radiance  observations  obtained  by  a  brassboard 
version  of  the  Infrared  Temperature  Profile  Radiometer  (ITPR)  from  the  NASA 
Convair  990  aircraft  during  June  1970.   Results  presented  were  obtained  from 
observations  during  a  flight  over  the  Pacific  Ocean  (from  Fairbanks,  Alaska 
to  San  Francisco,  California)  in  which  a  cold  front  was  encountered  near 
weather  ship  "Papa"  (50°N,  145°W).   These  results  demonstrate  that  spectral 
radiometric  measurements  are  very  useful  for  obtaining  a  detailed  depiction 
of  the  atmosphere's  meteorological  state  from  an  aircraft  observation 
platform . 


"Presented  at  the  Conference  on  Atmospheric  Radiation, 
Ft.  Collins,  Colorado,  August  7-9,  1972. 


2.   THE  AIRBORNE  ITPR 


The  1970  airborne  ITPR  experiment  has  already  been  discussed  in  detail  by 
Smith  et  al.,  (1972).   Briefly,  the  ITPR  measured  the  upwelling  radiance  from 
the  underlying  earth  and  atmosphere  in  five  different  spectral  regions;   one 
in  the  rotational  water  vapor  band  near  19um,  three  in  the  15-um  CO2  band, 
and  one  in  the  atmospheric  window  at  Hum.   The  detailed  spectral  response 
characteristics  of  the  ITPR  channels  are  presented  in  figure  1. 
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-Measured  transmittances   of  the   spectral    filters   in 
ITPR  channels    1   through   5 


3.      RADIANCE   SPECIFICATION  OF  ATMOSPHERIC  VARIABLES 

A.      Atmospheric  Temperature  Profiles 

The  radiance  measured  by  the   ITPR,   Iv  (pa),   is   related  to  the  atmospheric 
temperature  profile  through  the  radiative  transfer  equation 


MPa)=IVPo)Tv(Pa 


/Po 
Pa 


Bv(p)      dTv(Pa'P) 
dp 


dp. 


(1) 


B  (p)  is  the  Planck  radiance  given  by 

Bv(p)=ClV3/  rexp(C2v/T(p))-ll   ,  (2) 

where  C^  and  C2  are  constants  of  the  Planck  function,  and  T(p)  is  the 
temperature  at  pressure  p.  tv   (p  ,p)  is  the  transmittance  of  the  atmosphere 
between  the  pressure  p  and  the  aircraft  pressure  level  pa,  and  pQ  denotes  the 
surface  pressure. 

The  solution  of  (1)  for  T(p)  has  been  discussed  by  numerous  authors  and 
most  of  the  viable  methods  have  been  recently  summarized  by  Fleming  and  Snrth 
(1972).   For  the  aircraft  inference  problem  the  "statistical  regression" 
solution  has  been  adopted.   In  this  solution 

n 

Br(p)=Br(p)+  J2    ai<Pa>P)   [Ri(Pa)-Ri(Pa>]   >  <3> 

i=l 

where  Br(p)  is  the  Planck  radiance  corresponding  to  a  reference  frequency  r, 
chosen  to  be  intermediate  to  the  observation  frequency  range,  and  R^(pa)  is 
the  measured  radiance  for  the  ith  frequency  1^ ,  normalized  to  the  reference 
frequency  by  using  the  relation 

Ri  =  Clr3  '    ["exP(C2r  /   T?)"1]  >  (4) 

where  the  brightness  temperature  T?  is  given  by  the  inverse  Planck  function 


T?  =  C2  v£/  In   ["(^  Vi3/  I.)  +1 


(5) 


The  barred  quantities  denote  the  means  of  the  statistical  sample  used  to  de- 
rive the  regression  coefficients   ai(pa,p).   Once  Bp(p)  is  determined,  the 
temperature  profile  is  calculated  using  the  inverse  Planck  function  in  the 
form 


T(p)  =  C2r/  In  j  (Cjr3/   Br(p))  +  ll 


(6) 


In  this  study  the  regression  coefficients  were  determined  from  radiances  cal- 
culated for  the  heterogeneous  sample  of  106  atmospheres  given  by  Wark  et  al. 
(1962).   The  transmittance  functions  needed  to  calculate  the  ITPR  radiances 
were  taken  from  Smith  et  al.  (1972).   Random  errors  on  the  order  of  0.25% 
were  added  to  the  calculated  radiances  to  simulate  actual  radiance  observations. 
The  regression  coefficients  were  obtained  for  aircraft  pressure-altitudes  of 
950  mb  to  200  mb  in  increments  of  50  mb.   For  each  aircraft  pressure-altitude, 
regression  coefficients  relating  the  aircraft-level  temperatures  and  the 
"observed"  radiances  to  the  temperatures  of  various  levels  were  determined, 
using  the  standard  least-squares  solution. 

Figure  2  shows  the  resulting  standard  error  of  estimate  associated  with  the 
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Figure  2. --Standard  error  of  estimate  of  radiance  derived  temperature 
profiles  as  a  function  of  aircraft  flight  level 

vertical  temperature  profiles  for  the  dependent  sample  as  a  function  of  the 
aircraft  flight  level.   As  expected,  the  standard  error  of  estimate  for  the 
temperature  of  a  particular  pressure  level  beneath  the  aircraft  decreases 
with  decreasing  aircraft  altitude  (increasing  pressure-altitude),  especially 
for  the  pressure  levels  near  the  aircraft  pressure-altitude.   The  standard 
errors  of  estimate  are  largest  in  the  lower  troposphere  (800-950  mb)  where 
surface  inversions  occur,  and  above  300  mb  where  the  tropopause  usually 
exists.   The  temperature  errors  are  a  minimum  at  the  aircraft  flight  level 
because  the  measurement  at  this  level  is  direct,  and  are  also  minimal  at  the 
earth's  surface  because  the  window  channel  provides  a  nearly  unique  observa- 
tion of  the  surface  temperature. 

B.   Total  Precipitable  Water 


To  estimate  the  total  precipitable  water  vapor,  UQ ,  from  the  19- urn  water 


vapor  channel  radiance,  I,  the  standard  iterative  equation 


UJ+1=  rji   + 
o     o 


Uj  -  U  j"1 
o    o 

ft  _  £  j-i 


(I 


$ 


(7) 


A 


is  used.   In  (7)  I  denotes  a  calculated  radiance  and  the  superscript  denotes 
the  iterative  step.   The  precipitable  water  vapor  profile  required  to 


calculate  the  radiance  at  each  Iterative  step  is  given  by 


o 
' '  o 

where  A  is  used  at  its  climatological  value  of  3.0  (Smith  1966). 

C.   Clear  Column  Radiance 

When  clouds  exist,  the  radiance  propagating  from  the  clear  air  columns 
must  be  specified  from  the  measured  radiance  distribution  before  soundings 
down  to  the  earth's  surface  can  be  calculated.   From  an  aircraft,  one  may 
measure  many  spatially  independent  radiances  over  a  short  distance  along  the 
flight  track.   Consequently,  the  variation  of  the  observed  radiance  over  this 
short  distance  will  be  due  mainly  to  variations  in  the  cloudiness  within  the 
instrument  field  of  view  rather  than  to  variations  in  atmospheric  temperature. 
If  the  cloud  height  is  constant  over  the  distance  traversed,  so  that  the 
radiance  variations  are  due  to  variations  in  cloud  amount,  the  spatial 
variation  of  the  radiance  in  any  spectral  channel  will  be  linear  with  respect 
to  the  variation  of  the  radiance  in  any  other  simultaneously  observed  spectral 
channel.   This  linear  function  can  be  used  to  specify  the  clear  column  radi- 
ances from  the  observed  cloud  contaminated  radiance  distribution  (Smith  1968, 
Luebbe  1971). 

To  show  this,  consider  the  following  equation  governing  the  measured 
radiance  I(v)  when  a  cloud  fills  the  fraction  N  of  the  instrument's  field  of 
view: 

I(v)  =  NIcd(v)+(l-N)Ic(v)  (9) 

Here  s  lc(j  (v)  is  the  total  radiance  propagating  from  all  of  the  cloudy  columns 
of  air  and  Ic  (v)  is  the  radiance  arising  from  all  of  the  cloudless  columns 
of  air  within  the  instrument   field  of  view.   Differentiating  (9)  with 
respect  to  the  spatial  coordinate,  S,  assuming  IC(j  (v)  and  Ic  (v)  are 
constant  yields 

dl(v)    f  "I  dN 


_=Llcd(v)-Ic(v)j_   .  (10) 

Thus  it  follows  that  the  relation  between  the  radiances  observed  simultaneous- 
ly, at  two  different  spectral  frequencies,  v^  and  Vj ,  is 

A^Vj.Vj)  =  Wvi>  -  MV   =   dI^    ,  (ID 

I  ,(v.)  -  I  (v.)     dl(v.) 

cd  ]     c  ]         j 

where  A,  (v., v.)  is  a  constant.   Integrating  (11)  gives 

Kv.)  =  AQ(vi,vj)  +  A1(v1,vj)  I(v.)  (12) 


The  coefficients  A  and  A,  can  be  determined  from  a  set  of  simultaneous 
observations  by  standard  least-squares  regression  techniques.   Then,  given 
clear  column  radiance  for  any  spectral  channel  v-   the  clear  column  radiance 
for  any  other  spectral  channel  v^  can  be  predicted  from  (12). 

In  practice,  the  clear  column  radiance  for  the  window  channel  is  used  to 
predict  the  clear  column  radiance  for  the  water  vapor  and  CO2  channels  (see 
fig.  3).   The  clear  column  radiance  for  the  window  channel  can  be  obtained 
from  statistical  histogram  analyses  of  the  window  radiance  observations 
(Smith,  Rao,  Koffler  and  Curtis  1970)  if  a  sufficient  number  of  cloudless 
fields  of  view  exist.   If  no  completely  cloud-free  fields  of  view  are  avail- 
able, the  clear  column  radiance  can  be  specified  by  horizontal  interpolation 
or  from  an  estimate  of  the  surface  temperature,  i.e. 9  Ic  =  B  (T  ~  )  +  small 
atmospheric  correction. 

Figure  4  shows  the  mean  and  the  range  (maximum  and  minimum)  of  radiances 
observed  from  about  35,000  feet  within  50  n.mi.  sections  of  an  Alaska-to- 
California  flight.   Each  section  consists  of  100  spatially  independent 
observations.   The  clear  column  radiance  profile  for  the  window  channel 
(899  cm--'-)  was  obtained  by  spatial  interpolation  of  the  window  radiances 
obtained  from  cloud-free  fields  of  view  because  the  horizontal  variation  is 
small  over  the  sea.   The  clear  column  radiances  for  the  water  vapor  and 
carbon  dioxide  absorption  channels  were  then  calculated  for  each  leg  using 
the  regression  procedure  outlined  above.   These  clear  column  radiances  were 
then  used  to  calculate  the  temperature  and  water  vapor  profiles  (shown  later) 
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Figure  4.--ITPR  measured  radiances  and  calculated  clear  column  radiances 

from  flight  over  Pacific  Ocean,  June  12,  1970.   The  heavy  dashed 
line  represents  interpolated  values  between  data  at  first  point 
and  at  2030  Z. 

D.   Cloud  Height 


Assuming  zero  cloud  reflectivity  for  the  frequency  v.,  it  can  be  shown  that 


I  ,(v)  = 

cd 


(v)  iBcd(v)  +  ]_i-  e  (v) 


I  (v) 


(13) 


where  E     (v)  is  the  emissivity  of  the  cloud.   The  "black"  cloud  radiance  is 


IBcd(v)=Bv(pc)Tv(p,spJ   -  f   Bv(p)  dVPa'p) 


a  c 


/ 


dp 


(14) 


dp 


where  p   is  the  effective  cloud  pressure.   Substituting  (13)  into  (11)  gives 


dl(v±) 


dl(v  ) 


(v.) 


(V) 


IR  ,(v.)-I  (v.) 

Bed  i   c  i 

U   ,(v  )-l  (v  ) 

Bed   ]    c   ] 


(15) 


If  we  choose  two  frequencies  which  have  approximately  the  same  cloud 
emissivity  (e.g.,  channels  in  the  wing  of  the  15-ym  C09  band)  then  it  follows 
that 


A1(vi,v ,,pc)  =   IBcd(vi-Pc)  -  Ic(vi) 


I   ,(v.,p  )  -  I  (v.) 

Bed   J   c     c   3 


(16) 


Given  the ^temperature  profile,  as  calculated  from  the  clear  column  radiances, 

I   ,  and  I   can  be  calculated  using  (1M-)  as  a  function  of  cloud  pressure. 
Bed      c  to  v 

Hence,  using  (16),  the  radiance  slope  Aj_  (pc)  can  be  calculated  as  a  function 

of  cloud  pressure.   Then  the  actual  cloud  pressure  can  be  specified  from  the 

observed  slope  A, (v., v.)  using  the  iterative  solution 


^+1   3  pj 
c     ^c 
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Ai(pc)-Al(Pc_1) 


A  (v.,V.)-A  (pD 
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(17) 


which  converges  for  Aj_(v£,v-:)  f   0.   Also,  the  pressure  of  a  cloud  within  the 
field  of  view  of  a  single  observation,  within  the  sample  used  to  determine 
the  clear  column  radiances,  can  be  specified  by  letting 


I(vi}  -  Ic(vj) 

vvy =  Kv.) .  1  (v.) 


(18) 


E.  Effective  Cloud  Amount 


The  effective  amount  of  cloud,  N   (equal  to  N£(v)),  within  the  instrument 
field  of  view  can  then  be  specified  from  one  of  the  two  channels  used  to 
calculate  the  cloud  pressure  using  the  equation 


N"  =  N£(v)  = 


I(v)-Ic(v) 


(19) 


where  I  (v)  is  either  the  mean  of  a  sample  of  observations  or  an  individual 
observation,  depending  upon  which  was  used  to  specify  pc  and  consequently 

^cd^)- 


4, 


RESULTS 


Figure  5  shows  a  cross-section  of  atmospheric  temperature  obtained  while 
flying  at  33,000  feet  over  a  Pacific  cold  front.  The  temperature  profiles 
were  obtained  from  the  clear  column  radiances  shown  in  figure  4,  using  the 
statistical  regression  solution  presented  earlier.  The  isolines  represent 
departures  from  the  level-mean-values.  As  shown,  the  variation  of  temper- 
ature across  the  front  could  be  diagnosed  in  this  case  because  there  were 
sufficient  breaks  in  the  cloudiness  associated  with  the  front. 


Figure  6  compares  the  1000-  to  500-mb  and  the  500-  to  250-mb  thickness 
values  obtained  from  the  infrared  derived  temperature  profiles  with  thickness 
values  derived  from  simultaneous  microwave  measurements  (NEMS)  by  Rosenkranz 
et  al.  (1971)  and  observed  by  the  Weather  Ship  "4YP"  (Papa)  radiosonde.   As 
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Cross  Section  of  Aircraft  ITPR  Deduced  Temperature 
over  Pacific  Cold  Front,  June  12,  1970 


2200Z 


53N      53N      52N       51 N       51 N       51 N       50N      SON      50N      49N      49N 
I50W    150W    150W    149W     149W    U8W     147W    145W    U4W    143W    142W 

Flight  Track 


Figure  5. — Cross  section  of  ITPR  deduced  temperature  profiles  over  Pacific 

cold  front.   Isolines  show  the  departure  from  the  level-mean-values. 
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Figure  6. — Thickness  values  (1000-  to  500-mb  and  500-mb) ,  derived  from  infra- 
red measurements,  microwave  measurements,  and  4YP  radiosonde  data, 
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shown,  there  is  relatively  good  agreement  between  the  infrared  and  the  micro- 
wave data,  although  both  display  some  minor  discrepancies  with  the  conventional 
radiosonde  data  at  the  location  of  the  weather  ship. 

Figure  7  compares  the  total  precipitable  water  vapor  derived  from  the  infra- 
red water  vapor  channel  radiance  data  with  that  obtained  from  microwave  and 
radiosonde  data.   Here  we  see  large  discrepancies  between  the  total  precipit- 
able water  vapor  estimates.   The  systematic  discrepancies  between  the  two 
estimates  are  clearly  related  to  the  temperature  discrepancies  shown,  in 
figure  6.   These  result  because  an  over  estimate  (or  an  under  estimate)  of 
temperature  requires  an  over  estimate  (or  under  estimate)  of  water  vapor  to 
satisfy  a  given  water  vapor  channel  radiance  measurement.   The  differences  in 
the  relative  variations  through  the  frontal  zone  probably  are  due  to  the  fact 
that  the  infrared  estimate  is  only  for  the  clear  air  column  while  the  micro- 
wave estimate  includes  the  water  vapor  contained  within  the  clouds.   This  is 
verified  by  the  two  "Papa"  raob  values ,  one  of  which  was  observed  in  cloud- 
less air  while  the  other  was  observed  during  an  ascent  through  the  frontal 
clouds.   The  inability  of  an  infrared  sensor  to  sound  through  clouds  is  a 
severe  limitation  in  the  sensing  of  total  precipitable  water. 

Figure  8  shows  cloud  heights  and  amounts  estimated  using  two  CO2  channels 
(as  outlined  above)  and  those  determined  directly  from  window  radiance 
measurements  (i.e.,  the  pressure  at  which  the  air  temperature  was  equal  to  the 
window  channel  brightness  temperature),  and  the  comparison  of  both  with  visual 
observations.   The  comparison  of  the  two  estimates  shows  that  the  window 
temperature  deduced  heights  are  consistently  lower,  and  exhibit  a  large  degree 
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Figure  7. — Total  precipitable  water  vapor  derived  from  infrared,  microwave, 
and  4YP  radiosonde  data. 
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ances  using  732  cm"-1-  and  750  cm--1-  C02  channel  observations  (solid 
curve),  "window"  channel  radiance  (solid  bar)  and  visual  observations 
(symbols).   Shown  below  is  the  estimated  cloud  amount  derived  from 
radiances  measured  in  the  732  cm-1  and  750  cm-1  channels. 

of  variability  even  in  totally  overcast  situations.   This  variability 
apparently  is  due  to  variations  in  the  opacity  of  the  cloud.   The  measurements 
in  the  CO2  band,  on  the  other  hand,  provide  qualitatively  good  estimates  of 
the  cloud  heights,  in  comparison  with  corresponding  visual  observations,  even 
when  the  field  of  view  is  not  completely  filled  with  cloud.   This  is  demon- 
strated by  the  estimates  obtained  near  the  beginning  of  this  flight  track.   In 
thick  overcast  regions ,  the  CO2  channel  estimates  display  agreement  with  the 
maximum  heights  obtained  from  window  radiances. 

5.   CONCLUSION 

This  study  shows  that  multi-spectral  infrared  radiometers  aboard  aircraft 
can  be  a  very  effective  meteorological  observing  tool.   We  propose  that  this 
observing  capability  be  used  to  augment  conventional  observations  for  future 
sub-synoptic  scale  meteorological  research  programs,  such  as  the  forthcoming 
GARP  Atlantic  Tropical  Experiment. 
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